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We consider that the pre- inflation era is radiation-dominated, transiting smoothly to the inflation- 
ary era. We work out in detail the dynamics of inflaton fluctuations across the phase transition and 
the proper choices of initial vacuum states. It is found that this phase transition can suppress long- 
wavelength quantum fluctuations of inflaton. This may attribute to the large-scale CMB anisotropy 
a lower power than predicted in the standard ACDM model. In constraining this transitional effect 
by WMAP anisotropy data, we use the WMAP best-fit scale-invariant ACDM model with the den- 
\ sity power spectrum replaced by the one found in this work. We find that the transition occurs at 

least about 10 e-folds before the comoving scales comparable to our present horizon size cross the 
1 Hubble radius during inflation. 
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PACS numbers: 98.80.Cq, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

The inflationary model can explain the observed spatially flat and homogenous Universe. Moreover, quantum 
fluctuations during inflation can give rise to primordial density fluctuations with an almost scale-invariant power 
spectrum, which is consistent with the recent WMAP data on cosmic microwave background anisotropies [l| . Among 
various cosmological models, the ACDM model is a quite good one - it can fit CMB data well. However, the CMB 
anisotropy measurements made by WMAP have shown that the amplitude of the quadrupole is lower than expected 
for the ACDM model, although the statistical significance of such an anomaly is not large Recently, there have 
been many proposed solutions for this anomaly, mostly based on some new ingredients in the generating process of 
density perturbation (2j. 

A general assumption of the ACDM model is that it considers only a vacuum-dominated inflationary epoch, and 
that it does not consider the possibility of a pre-inflation era. For example, a radiation-dominated era has taken place 
before inflation. In fact, this phase transition has been proposed in the context of spontaneous symmetry breaking 
to provide with a mechanism for inflation ||. In this paper, we will take into account the effect of a pre-inflation 
1 radiation-dominated era. As expected, if inflation lasts very long such that the present Universe only occupies a 
tiny portion of the inflated region, the effect of the phase transition will be negligible. However, if it happens that 
the present Universe is just comparable to the size of the inflated region, the pre-radiation era may leave imprints 
^7"; ' on large-scale structures of the Universe. Here we will study the effect of the phase transition to the large-scale 
J> , CMB anisotropy. Recently, this effect has been discussed by Powell and Kinney who considered large- wavelength 
inflaton fluctuations in the pre-inflation radiation-dominated era and found a suppression of power on large scales as 
a result of the choices of vacuum states. (See also Ref. [j| for a similar consideration.) In addition, Contaldi et al. 
in Ref. have found a similar effect caused by a pre-inflation epoch that is dominated by the kinetic energy of the 
inflation. However, the issues about the dynamics of inflaton fluctuations across the phase transition and the proper 
choices of initial vacuum states have not been fully investigated. Here, we will study these issues in detail. 

Let us simply consider that the Universe has two components during the phase transition - one is radiation and the 
other is vacuum energy. We assume that the pre-inflationary Universe is in a radiation-dominated phase and then the 
Universe evolves to a vacuum-energy-dominated or de Sitter phase. Although the solution for this transition has been 
obtained for many years [f| 0, B @| ■> there is no attempt in varying the amounts of radiation and vacuum energies of 
the Universe to see the effect on CMB anisotropy. This serves a way to test the existence of an unknown pre-inflation 
phase. Also, it provides a constraint on the period of inflation which needs to be consistent with the observed CMB 
anisotropy spectrum under such a pre-inflation condition. 

In this paper, we use Vilenkin and Ford's transition solution [6| to compute the power spectrum of quantum fluc- 
tuations of the inflaton field. Since the pre-inflation phase is radiation-dominated, we use the radiation-dominated 
solution to determine the boundary conditions. This allows us to set the form of the seed for the quantum fluctuations 
before the de Sitter phase begins. By introducing this pre-inflation initial condition, we can calculate the power spec- 
trum of the horizon-crossing modes near the transition region. Then we input the power spectrum to the CMBFAST 
numerical code [l(| to generate the CMB anisotropy power spectrum. We vary the amounts of radiation and vacuum 
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energies and tune the duration of inflation to obtain different results. In order to compare the results with the ACDM 
model, we do the chi-square fitting of the generated CMB anisotropy power spectra to the WMAP data and calculate 
the chi-square values. 

In next section, we will derive the mode equation of the quantum fluctuation during the phase transition. The 
numerical results for the mode equation and the respective CMB anisotropy power spectra will be given in Sec. IIIII 
Sec. IIVI is our conclusion. 



II. THE MODEL 



We consider a flat Robertson- Walker metric, with a scale factor a(t) : 

ds 2 = g^dy^dy" = dt 2 - a 2 (t)dy 2 . 
The total energy-momentum tensor is given by 

V = + Vb<V , 

where T^™^ is the trace-free energy momentum tensor for radiation and Vq is constant vacuum energy. Hence, the 
Einstein evolution equation is 



(1) 



(2) 



H 2 



1 da 
a dt 



3M 2 pl ' 



where M Pi = (8ttG)-^ 2 = 2.436 x 10 18 GeV is the reduced Planck mass and p is given by 



(3) 



(4) 



with Vq = 3BMq of order of the grand unification energy scale. A and B are dimensionless constant parameters. 
Defining A 2 = M%jM pl , the solution of Eqs. © and g]) is found as [1] 



a(t) = 



A\ 3 r 

B 



sinh [2y f BAt 



(5) 



The asymptotic form in the radiation-dominated phase is 



At < : a(t) - y/2Ai(At)i 



(6) 



The asymptotic form in the inflationary phase is 



At > -B-* , a(t) - T 



B 



,v / SAt 



(7) 



To study the transition from the pre- inflation radiation-dominated phase to inflation, we put Eq.© into the 
equation of motion for the inflaton quantum fluctuation. To a good approximation, we consider the zero-mass case. 
The equation of motion for the quantum fluctuation of inflaton is given by 



a<fi(t,y) =0, 



(8) 



where 



□ = — 



-ft 



The Fourier component of tp(t, y) is given by 



•dy. 



(9) 



(10) 
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Put Eq. (JTOj) into Eq. ®, then Eq. © becomes 

Mt) + 3-Mt) + (K) = o , kit) = . (ii) 

To solve Eq. (fTTj) with a(i) given by Eq. ([5]), we change the derivative of an independent variable from the time t 
to the scale factor a(t). As we want to study the effect of the transition around the critical point a c , where a c is given 
by Aja\ — B, we rewrite Eq. (fTTj) in terms of x = a — a c . Then, Eq. (fTTj) becomes 



[B(x + a c ) 4 + A] <fi(x) 



AB(x + a c ) 3 + 



2A 



x + a c 



<j>' k (x)+k 2 4> k {x) = Q > ti(.x) 



_ rf<fa (x) 



(12) 



Note that when a(t) is small, the Universe is in radiation-dominated era. Therefore, we choose the boundary conditions 
for Eq. (fUl as 



.-(*) 



1 1 



a V2fc 



2ka 2 



2Aa 



a ika/ y/~A 



(13) 
(14) 



Solving Eq. (|12[) with above boundary conditions gives the solution for a particular k mode. The relation k = aH 
determines the value of k for the mode that crosses the horizon at a = at, where H is the Hubble parameter given 
by Eq. ©. This gives from Eq. © that 



A: 
A 



"2 + Ba % 



(15) 



Since we are concerned with the fc-modes that leave the horizon shortly after inflation begins, we will consider only 
the solutions with k > k c , where k c — K^/2Ba c . The power spectrum of the quantum fluctuations of the inflation is 
then given by 



P fc 1/2 cx fe 3/2 </> fe (z fe ) , 



(16) 



1/2 

where Xk = au - a c - In the following section, we compute power spectra P k produced when and after inflation 
begins (x > 0) by numerically solving Eq. (TT2|) for different values of A and B. Then, we put these results into the 
CMBFAST code [l(| to generate the CMB anisotropy power spectra. 



III. NUMERICAL RESULTS 



So far, we have not specified the values of the parameters A and B. Before doing the numerical calculation, let us 
assess the possible range of the values of A and B. 

From Eq. the vacuum energy that drives inflation is given by 

V a = MI<kB . (17) 

Here we consider Vq as the potential energy of an inflaton potential V(ip) during inflation. The measurements of the 
CMB anisotropy made by WMAP has put a constraint on V{ip) given by [Il| 

V* = 0.0265e3 M PU (18) 

where the slow-roll parameter 

e = {M 2 Pl /2){V' /V) 2 < 0.033 (19) 
at 95% confidence level. Therefore, we have Vq ~ V and thus 

< 0.0086M Pi /M G . (20) 
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If we take M G = 2.1 x 10 16 GeV, then B < 1. 

In the following, we choose B = 1 and B = 0.1 to calculate the power spectra of inflaton fluctuations P k in Eq. (|16[) . 
Although we have introduced the slow-roll parameter e, we will assume a pure de Sitter space during inflation and 
neglect the small slow-roll corrections to the density power spectrum. This would not affect our interest here since 
the effect of a pre-inflation radiation-dominated epoch is mainly on large angular scales. Then, we use the resulting 
density power spectrum and the other cosmological parameters in the WMAP best-fit scale-invariant ACDM model 
(i.e. the density power spectrum has a scalar spectral index n s = 1) |l( to generate the CMB anisotropy power 
spectra. Note that while the observation gives the constraint on B, the value of A is undetermined. In fact, different 
values of A introduce different boundary conditions as shown in Eqs. (|13p and (|14p . Therefore, we pick a certain value 
of A to see how the value affects the CMB anisotropy power spectrum and study under what conditions the results 
are compatible with present CMB observations. In order to give a quantitative comparison with the CMB data, we 
use the chi-square fitting, x 2 = ^ z — T^)C~^(D^ — T^), and the chi-square value as a measure, where D\ is the 
measured zth band-power, is the theoretical value, and C ltJ is the width of the error bar in measurement. 

We present the following two cases: 

(1) B = 1 with A = 7, 5, 1, 0.1. In this case, the transition takes place at a c = 1.63 — 0.56 with k c /A = 2.30 — 0.80. 

1 /2 

Figure Ufa) shows the power spectrum of inflaton fluctuations P k for each set of B and A. It is apparent that 
the large-scale or small-/c power is suppressed by the presence of the pre-inflation radiation-dominated epoch. This 
suppression is similar to that found in Ref. However, there are two main differences. For large-fc modes, we have 
assumed a standard inflationary scale-invariant power spectrum with n s = 1 while they have taken the spectrum 
with n s — 0.951. For small-k modes, our spectrum increases monotonically with k whereas their spectrum undergoes 
oscillations; however, both spectra rapidly approach the inflationary spectrum. This oscillatory behavior, also found 
by Contaldi et al. in Ref. [2J (the result from an analytic computation as shown in Fig. 1 of their paper) and Enqvist et 
al. in Ref. @, is actually an artifact resulted from the discontinuity of Ricci curvature at the abrupt phase transition. 
This artificial particle production can be avoided by considering for example a C°° function for the scale factor 
a{t) [12|, as we have worked with in Eq.([5]). 

Besides, we must be very careful about choosing the initial point while doing the numerical calculation. For 
example, in our case with B = A = 1, the phase transition point is at a c = 1. We have obtained a smooth power 
spectrum (see Fig.QJa)) by choosing the initial a, to be ai — 0.0001 (which means that it is in the very early Universe; 
therefore, it is purely radiation-dominated and also very far from the phase transition point a c = 1). But, if we tune 
the initial <Zj to be close to the phase transition point, say dj = 0.3 , the oscillation begins to show up. And if we 
tune the initial a, closer to the phase transition point, for example, a, = 0.6 , the oscillation becomes larger than 
that in the dj = 0.3 case, as shown in Fig. [Kb). This test explains that if the initial <Zi is not in the proper region 
(proper region means the radiation-dominated region and that is the right region for the radiation-dominated solution 
in Eq. (|13p . while improper region means the phase transition region which is not a radiation-dominated region and 
therefore the solution for <j>k(x) is not Eq. (|13|) anymore). This improper choice for the initial a% indeed leads to an 
oscillating power spectrum. However, if the initial a, is carefully chosen in the proper region, it does not lead to an 
oscillating power spectrum. Note that here we have used the scale factor a as the evolution variable. If we have used 
the conformal time, we would have obtained an oscillating power spectrum which approaches the flat inflationary 
spectrum at large k. This explains why the power spectrum obtained by Contaldi et al. in Ref. [2| (see Fig. 2 of their 
paper) is oscillating although it is obtained from an exact numerical evolution. In their paper, they have indeed used 
an improper choice of the initial condition which is too close to the phase transition. (They have chosen an initial 
kinetic energy of the inflaton only about 100 times of the potential energy. Since the kinetic energy of the inflaton 
decreases as a -6 , both energies rapidly become comparable. To obtain a correct power spectrum, they should have 
used a much larger initial kinetic energy.) 

1 12 

The power suppression would affect the CMB anisotropy power spectrum. For a spectrum P k , we choose a value 

of z = fco.os/A corresponding to the physical scale of 0.05Mpc _1 and then input this P]J 2 into the CMBFAST code 
to generate the CMB anisotropy power spectrum C;. This value of z also corresponds to an exponential expansion 
with the number of e-folds from the start of inflation given by 

N z = \n(a k /a c ) ~ lnz. (21) 

The results are shown in Fig. [51 Since the measured low-/ multipolcs of the CMB anisotropy power spectrum are 
slightly different in the three-year WMAP data (WMAP3) [l]| from the one-year data (WMAP1) Q, we compare 
the chi-square fitting by using WMAP3 data with that with WMAP1 data. The results are summarized in Tables [J 
and OH 

The x 2 value of the ACDM model with a scale-invariant power spectrum using WMAP3 data is 47.09. From 
Table [H we can see that the x 2 value for the model decreases to about 48 as z increases to 3 x 10 4 , corresponding 
to N z ~ 10 e-folds after the start of inflation. This can be reflected, for example in Fig. ^Bp, by those five curves for 
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z = 120, 300, 500, 1000, 3000 rising in turn from down to top. As we further increase the value of z, the resulting Ci 
approaches to that of the ACDM model with a scale-invariant power spectrum. Thus, it can be concluded that the 
WMAP3 data does not prefer the presence of a pre-inflation radiation-dominated epoch, although it is not excluded. 
If we assume that the comoving scales comparable to our present horizon size cross the Hubble radius during inflation 
when there are N cm f, ~ 60 e- folds from the end of inflation [la ] , then the radiation-dominated epoch would take place 
at least N = N z + N cm b ~ 70 e-folds from the end of inflation. 

The case is different when using WMAP1 data since the measured one-year low-/ multipoles are slightly smaller than 
the three-year results. The x 2 value of the ACDM model with a scale-invariant power spectrum using WMAP1 data 
is 64.99. From Table HT1 we can see that the x 2 values for all other models are slightly lower than this value at z ~ 10 4 , 
corresponding to N z ~ 9 e-folds after the start of inflation. Therefore, the WMAP1 data would, albeit a rather low 
significance, indicate that the universe would be actually predominated by radiation when N — N z + N cm b ~ 69 
c-folds from the end of inflation. 

(2) B = 0.1 with A = 7,5, 1,0.1. The transition takes place at a c = 2.89 - 1.00 with k c /A = 1.29 - 0.45. We do 
not show the power spectrum of inflaton fluctuations Pk for this case which is indeed quite similar to the previous 
case. Figure [3] shows the Ci power spectra and their x 2 values are given in Tables IIIII and IIVI In Table IIII1 the 
X 2 value decreases rapidly to about 47 at z ~ 3000; however, at large z all models approach to the ACDM model 
with a scale-invariant power spectrum. In Table ITVl the % 2 values for all the models reach a minimum at z ~ 1000, 
which is lower than the x 2 value for the ACDM model with a scale-invariant power spectrum. This corresponds to 
an inflationary expansion with e-folds of N z ~ 7 after the start of inflation. Because the x 2 values decrease with z in 
the B = 0.1 case faster than the B = 1 case, the anisotropy power spectra are closer to each other and the B = 0.1 
case reaches the standard ACDM model at a smaller z value. It means that inflation in the B = 1 case lasts longer 
than the B = 0.1 case in order to be consistent with the CMB anisotropy power spectrum observed today. We have 
assumed that inflation is driven by a nearly flat potential and thus neglected any spectral index. As we can see in 
Fig. [3] there is a certain degeneracy in the CMB power spectrum between the variation of the spectral index and the 
existence of a pre-inflation radiation-dominated era, so a joint likelihood fitting of relevant cosmological parameters 
to WMAP data should be carried out before drawing any conclusion about pre-inflation physics. 

IV. CONCLUSION 

We have shown that the existence of a pre-inflation radiation-dominated phase can affect the CMB anisotropy 
power spectrum. This transition can give a low quadrupole value by adjusting the values for A, B, and z. The smaller 
is the value of z, the closer is the comoving scales of the present Universe to the transition when they cross out the 
horizon during inflation. This manifests as rising up or down for those low-/ C; of the z = 120 spectra in Figs. [Hand 
[3J Meanwhile, the high-/ C;, for example, in Fig. is significantly lower than the observed values. So, it is unlikely 
that the inflationary era lasts for only 60 e-folds. As expected, the effect of the radiation phase becomes weaker as 
the value of z increases. From the three-year WMAP data, using the WMAP best-fit ACDM model parameters, we 
have found that the phase transition occurs at least 70 e-folds from the end of inflation. Interestingly, the one-year 
WMAP data indeed prefers very slightly the pre-inflation radiation-dominated phase taking place as early as about 
70 e-folds from the end of inflation. Future CMB observations will provide more severe constraints on this transitional 
effect, and a joint likelihood fitting to CMB data for studying the degeneracy between the addition of a pre-inflation 
radiation-dominated era and other relevant cosmological parameters will certainly be necessary. 
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z = 120 


z = 300 


z = 500 


z = 10 3 


z = 3 x 10 3 


z = 10 4 


z = 3 x 10 4 


X 2 (B1A7) 


997 


367.2 


211.3 


103 


57.91 


49.27 


47.70 


X 2 (B1A5) 


1299 


443 


217 


102 


57.43 


49.16 


47.67 


X 2 (B1A1) 


3548 


449.4 


209.6 


98.37 


56.12 


48.87 


47.59 


X 2 (B1A0.1) 


1574 


411.6 


208.8 


99 


56.54 


48.96 


47.61 



TABLE I: Using WMAP3 data to the chi-square fitting of the CMB anisotropy power spectra for B — 1 and A — 7, 5, 1,0.1 
with different z values. Note that the chi-square value of the ACDM model with a scale-invariant power spectrum is 47.09. 





z = 120 


z = 300 


z = 500 


2 = 10 3 


z = 3 x 10 3 


z = 10 4 


z = 3 x 10 4 


X 2 (B1A7) 


639.50 


236.37 


146 


85.45 


65.93 


64.47 


64.72 


X 2 (B1A5) 


822.23 


280.6 


147.16 


84.83 


65.81 


64.48 


64.74 


X 2 (B1A1) 


2238 


283.35 


142.47 


82.90 


65.44 


64.50 


64.77 


X 2 (B1A0.1) 


992 


261.4 


142.2 


83.51 


65.55 


64.49 


64.75 



TABLE II: The same as in Table Q] but using WMAP1 data. Note that the chi-square value of the ACDM model with a 
scale-invariant power spectrum is 64.99. 





z = 120 


z = 300 


z = 500 


z = 10 3 


z = 3 x 10 3 


z = 10 4 


z = 10 S 


X 2 (B0.1A7) 


59 


56.88 


52.82 


49.83 


47.97 


47.43 


47 .11 


X 2 (B0AA5) 


67.3 


56.9 


52.59 


49.66 


47.89 


47.32 


47 .11 


X 2 (B0AA1) 


80.08 


56.19 


51.78 


49.08 


47.67 


47.23 


47.09 


X 2 {B0.1A0A) 


77.98 


55.18 


51 


48.63 


47.49 


47.2 


47.08 



TABLE III: x 2 values for B = 0.1 and A = 7,5, 1,0.1, using WMAP3 data. 
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z = 120 


z = 300 


z = 500 


z = 10 3 


z = 3 x 10 3 


z = 10 4 


z = 10 5 


X 2 (B0.1A7) 


67.13 


65.78 


64.8 


64.49 


64.66 


64.87 


64.99 


X 2 (B0.1A5) 


69.68 


65.75 


64.76 


64.49 


64.68 


64.87 


64.98 


X 2 (B0.1A1) 


74.51 


65.49 


64.6 


64.50 


64.75 


64.91 


64.99 


X 2 (B0.1A0.1) 


73.54 


65.21 


64.5 


64.51 


64.80 


64.95 


64.99 



TABLE IV: x 2 values for B = 0.1 and A = 7,5, 1, 0.1, using WMAP1 data. 
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FIG. 1 (b) 
i 1 1 1 1 1 r 




k/A 



FIG. 1: (a) Normalized power spectra of inflaton fluctuations with A = 7, 5, 1, 0.1 for a fixed B = 1. Inserted panel shows the 
spectra for small k. All curves start with k = k c corresponding to the fc-mode that leaves the horizon at the start of inflation, 
(b) Normalized power spectra of inflaton fluctuations with initial a,i — 0.001,0.3,0.6 for the case with B = A = 1. The phase 
transition point in this case is o c = 1. 
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FIG. 2: CMB anisotropy power spectra for B — 1 and A = 7, 5, 1,0.1. In each panel, the curves are induced by the with 
2 = fco.os/A = 120,300,500, 1000,3000 corresponding to 0.05Mpc~\ The z = 3000 curve almost overlaps with that of the 
ACDM model with a scale-invariant power spectrum. We normalize all the anisotropy spectra at the first Doppler peak. Also 
shown are the three-year WMAP data including error bars. 




FIG. 3 (c) B0.1A1 



FIG. 3 (d) BO.lAO.l 




FIG. 3: The same as in Fig. [2]but for B = 0.1 and A = 7, 5, 1, 0.1. 



